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x*his  report  corrects  errors  in  a previous  report  on  the  same 
subject  and  presents  a listing  of  a revised  and  improved  digital 
computer  program  for  finding  stress  distribution  in  a thin 
rotating  disk  with  nonuniform  heating.. 
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in 


A Computer  Subroutine  for  Stress 
Analysis  of  Rotating  Clous  - II. 

by  John  E.  Rrock 

Cased  upon  theory  developed  by  the  writer,  R.  E.  drown  developed 
a successful  computer  prolan  for  analysis  of  radial  and  circumfer- 
ential stresses  In  rotating  axisynr.etric  disks  of  variable  thickness 
having  an  axisynmetrical  thermal  strain  field.  The  writer  revised 
Brown* 3 program  so  as  to  invoke  a group  of  ancillary  subroutines 


which  luve  been  found  useful  in  another 
txowever,  much  unnecessary  and  confusing 
In  particular,  one  of  the  normalisations 
fail  in  the  quite  corncn  case  of  a disk 


application.  In  doing  so, 
normalisation  was  introduced. 

would  cause  the  analysis  to 
with  no  radial  loading  at  its 


outer  bo  uniary . All  tills  material  appears  as  Reference  1,  hereof. 

Referees  evaluating  a paper  based  upon  Reference  1,  called 
attention  to  these  faults  so  tint  the  program  has  been  rewritten. 

A listing  of  the  main  subroutine,  RODISK,  as  revise.!,  as  well  as 
listings  of  the  ancillary  subroutines  may  be  found  in  Appendix  A 
hereof,  rice  reader  will  note  that  other  changes  have  also  been  made 
resulting  in  somewhat  more  flexibility  of  application.  Employment  of 
the  revised  program  is  described  in  the  textual  material  which  appears 
at  the  beginning  of  the  listing. 

Appendix  9 contain  a revision  of  the  second  illustrative  example 
problem  of  Reference  1.  Thi3  problem  was  solved  for  various  values  of 
M ■ N-l,  the  number  of  equal  subdivisions  into  which  the  annular  radius 
b-a  is  divided  for  purposes  of  numerical  analysis  by  RODISK.  Also,  a 


number  of  different  values  of  KP(3)  were  used.  If  K?(3)  > 3,  it 3 
value  is  the  number  of  iterations  which  will  be  performed  by  R0DI5K. 
1*  KP( 3)  < 0,  iteration  will  continue  until  three  successive  values 
of  tne  unknown  parameter  3 determines  in  the  course  of  tne  analysis, 
satisfy  the  relation 


I Bi  — Bj  | + |92-B,|  + |3j-3i 
|B,j  + |b2|  + |b3' 


< 10 


,KP(3) 


.ve  also  determined  execution  time  by  use  of  the  library  subroutine 
IXCLOK,  executing  'under  CP-oms  on  the  I3M  360/67  at  the  W.  R.  Church 
Computer  Center  at  the  J&val  Postgraduate  School. 

We  found  that  execution  time  oer  iteration  is 


3 1.2  M + 5 (milliseconds) 


for  any  problem. 


Accuracy  was  evaluated  cy  dealing  with  problems  having  available 
analytic  solutions.  It  was  found  that  the  principal  limitation  on 
accuracy  is  determined  by  the  choice  of  subdivions,  the  integer  M * 

N-l,  so  that  there  is  a certain  inherent  error  regardless  of  how 
many  iterations  are  made.  This  error  depends  on  M,  of  course,  and  upon 
the  details  of  the  problem.  The  error  is  greatest  .near  the  inner 
radius  of  an  annular  disk,  and  is  large  if  the  ratio  a/b  is  snail. 
Fortuitously,  the  error  may  be  smaller  for  an  early  iteration  than 
for  a somewhat  later  iteration  but  this  is  not  practically  useful 
information.  For  the  problem  of  Appendix  3 hereof,  with  a/b  3 .165, 
we  find  the  results  given  in  Table  1,  (see  next  page). 

Thus,  for  example,  with  M • 20,  there  is  an  inherent  error  of 
about  IS  and  the  results  are  not  significantly  improved  by  iterating 


i 

I 


total 
time, 
secs . 

. 355 

• ^4- 

.32 
.90 

3.1 

Table  1.  Percent  error,  required  iterations,  and 
execution  time  for  problem  of  Appendix  3. 

more  than  eleven  times.  With  eleven  iterations,  the  solution 
is  returned  from  RODISK  in  0.32  seconds. 

The  significant  conclusion  is  that  the  execution  is  so  fast 
that  one  may  as  well  take  M * 100  (corresponding  to  M * 101,  the 
maximum  available  under  present  dimensioning)  and  iterate  many 
more  times  than  is  strictly  necessary.  'Taking  M • 101  and 
:<P(3)  * -3  gave  execution  in  3.7  seconds  with  31  iterations  and 
with  an  accuracy  of  0.004*  (In  the  oroblem  at  hand,  a_(a)  was 
soecifled  as  zero  and  the  program  gets  -1.14E-11  so  that  the  error 
here  is  "infinite'.'  Our  evaluation  of  0.0045  is  for  the  first 
position  rather  than  for  the  zeroeth. ) 

This  concludes  the  text  proper  of  the  present  recort.  however, 
we  take  advantage  of  thi3  opportunity  to  correct  errors  in  Reference 
1,  viz. : 

(1)  Page  3,  equation  12  snouid  read 


\1 

approx. 

limiting 

5 error 

approx 

iters . 

req'd. 

5 

16 

5 

10 

5 

7 

20 

1 

11 

40 

.1 

17 

100 

.01 

25 

m ■ */(n*-4vn+4)  ■ t/[(n-2)*+4(l-v)n] 
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Listing  of 


'urrout  ir.e 


30DISK 


mi  anciliarv 


subroutines 


SUBROUTINE  SOOI'K.  JOHN  E.  BROCK,  l MA  V 1978,  PJVISEO  1 AUGUST  1S7£.  RCDCCC1C 
THIS  IS  A SUBROUTINE  FCR  DETERMINING  SlOUt  AhO  £ IRCUMFERENT I AL  STRESS RC C CCC 00 
IN  AN  AXl  SYMMETRIC  THIN  ELASTIC  CIS*  HAVING  AN  AX  IS  VMM  ETR  IC  T HSR-*'  AL  ROOCCC2C 


IN  AN  AXI  SYMMETRIC  THIN  ELASTIC  0 1 S « HAVING  AN  AX  IS  VMM  ETR  IC  TMrRVAL 
STRAIN  FISLO  ANO  ROTATING  A 1 ANGULAR  VELOCITY  ONEGA  IRAQI ANE/SECCNC) 
ABOUT  THE  AXIS  OF  symmetry.  TwO  TYPES  OF  PFOBlEM  MAY  EE  TREATED: 


ABOUT  THE  AXIS  OF  SYMMETRY.  TwO  TYPES  OF  PPO&LE.M  MAY  £E  TRSATSO: 

TYPE  1:  ANNULAR  DISn  CF  IN'ICE  RACIUS  ARAC  ANO  OUTSICE  RACIUS 

efiAO.  the  .racial  stress  is  era  at  the  inner  raoils  anc 
srb  at  the  cuter  ajoius.  TMc  insice  raciue  must  be 
GREATER  Than  CERO 

TRTFE  2:  set  IC  OISK  HAVING  RAO  IAL  STRESS  SRB  AT  Ol/ISIOE  RAOIUS  BRAC. 


THE  USER  MUST  PROVIDE  A MAIN  PROGRAM  YhICh  CALLS  SLBROUTINE  ROCISK 
AFTER  IT  HAS  SLPPUSC  The  FOLLOWING  INFCRMAT  ION. 


C C * Mo 
RCCOOCSO 
RCCCCCtC 
RCCCCC7C 
RCOCCCEC 
RCCOOCSO 
RCOOO  10  D 

rcocc:  ic 

FCCCC  lie 


U)  N,  INTEGER.  IN-1)  IS  The  NUMBER  C F EGLAL  SUBDIVISIONS  INTO  WHICHRCC  CO  ISC 


THE  ANNULAR  RAOiuS  (BRAD  minuS  AR  A : J IS  0IV1CSC  FCR  CCMPU- 
TATICNAL  FuRPCSE  $•  THE  PRESENT  01  f>  ENS  ION  I NG  CAtf  ACCOMMODATE  N 
NOT  GREATER  THAN  101  . 

(2)  BRAD 

(3)  ARAC  (NOT  NECESSARY  FOR  PROBLEMS  OF  TYPE  2. > 

U>  SRB 

IS  J SRA  (NOT  NECESSARY  F0«  PROBLEMS  OF  TYPE  2.) 

(6)  POIS.  POISSON'S  RAT  10 

(7 I KP  (11*1,2*  INTEGER  TO  OEMOTE  PPCBlEM  OF  TYPE  1,2. 

(8)  KP(2>,  INTEGER  TO  PROVIDE  FCR  SKIPPING  WHILE  PRINTING 
OUTPUT.  FCR  EXAMPLE.  IF  N*1D1  AND  kF(2I*S,  GNLY  EVERY 
FIFTH  SeT  OF  VALUES  KILL  8=  PRINTEO:  1ST  t6TH, . . .,  96TH, 

ANO  10  1ST  . 

(9)  KP(3»,  INTEGER  SPECIFYING  NUHEER  OF  ITERATION  70  BE 


ROCCCl^O 
RCCCC1  IC 
ROOOO 1£J 
RCDCC170 
RCCCC16C 
ROD CC  ISC 
PC0CC2C3 


FQDC02  ID 

PCCCCi 2C 
RCDCCORC 
RCCOOIfD 
PCD  C C 160 
PCCCC2  IC 
SOCCC260 


«M3P  , ir  Uixu  u r iicnaiiui:  IU  oc 

PERFORMED.  USUALLY  kP(DI«10  IS  SUFFICIENT  FCR  ENGIN- 
EERING ACCURACY  . ALTERNATELY,  If  KPI3J  IS  A NEGATIVE 
INTEGER,  ITERATION  will  CONTINUE  UNTIL  THREE  SUCCESSI 
VALUES  OF  A PARAMETER  E.  DETERMINED  INTERNALLY,  ARE 
SUFFICIENTLY  CLOSE  Ai  COMPARED  TO  AN  EPSILON  ECUAL  TC 
TEN  RAISED  TO  THE  KF(2)  PCwER. 


R0CCC2SG 
PC0CC2CC 
RCCCC2  K 
FCCCC22C 
PCCCC22C 
RCCCC2AC 
RCCCC3  SC 
RCCCCStC 
RCCCC27C 
RCCCCStC 
RC0CC2SC 
RCDCOWCO 
PCDCCR  10 
RCDCCAiC 


(1X3)  KP(A).  IF  KP(4!*0,  CNLY  FINAL  ANSWERS  HILL  EE  PRINTED. 

IF  KP(4I*1,  A SEQUENCE  OF  ITERANT  VALUES  OF  ( hill  BE 
PRINTED  TO  l NCI  CATE  DEGREE  OF  CONVERGENCE.  IF  KPIA)>1. 
THERE  WILL  BE  NO  PRINTING  AT  ALL  WITHIN  ROOISK,  BUT  UPCN 
RETURN  KP(5I  HILL  CONTAIN  THE  NUMESR  OF  ITERATIONS  WM-ICH 
HERE  FERFCRMtO  SO  THAT  KP  ( 5 I MUST  BE  RESET  BEFORE  RCCISK 
IS  CALLEQ  AGAIN. 

(Ill  KP  ( 5 1 . WP(5»»0  CALJES  "ILNE  CUBIC  -SPLINE  INTEGRATION 

TO  BE  US E 0 • OTHERWISE  TRAPsiOICAL  INTECRATjCN  IS  USEC. 
<1.21  VECTCR  Xllwl,  v,*1.2,...,N,  CONTAINING  VALUES  C?  DISK 

THICKNESS  AT  eUUALLY  SPAcEO  ft  AC  I X FRCM  INSICg  TO  QUTSICE. 


VECTCR  X ( 2 , J I CONTAINS  VALUES  CF  Gamma  TIMS*  ONEGA 
S JUAR  EC  WhSSs  GAMMA  IS  (MASS  I CENS1TV  OF  Thi  MATERIAL. 

FOR  MOST  PROBLEMS  GAMMA  OCES  NOT  VARY  WITH  RACIUS  ANO 
ALU  ELEMENTS  OF  ThE  VECTOR  WILL  BE  THE  SAME. 

VECTOR  XI2.JI  CONTAINS  VALUES  CF  (JS  I ( ALPHA  I (TEE  I WHERE 
EE  IS  YOUNG'S  MG Cul US  , alpha  is  7 He  COEFFICIENT  OF  LINEAR 

thermal  expansion,  ano  tes  is  the  temperature  changc. 


RCrCCA2C 

PCDCCAAC 


RCCCOAIC 

RCOCCAtC 


RCCCCA7C 

PCOCCAtC 

PCCCCsSC 

RCCCC5CC 

RCCCCSTC 

RCCCCJCC 

RCCCCS2C 


Th®  MAIN  PRCGRAN  MUST  CONTAIN  THE  STATEMENTS: 
IMPLICIT  REAL  *8  <A-H,C-2I 


3 - 


I 


1 • 


*8*  fwSSS’“I"(,nUcrl“»« 


xi;o.  ion 

■R  \PL5l 


? 3¥t 

C T« 

5 &aYs  - 

c in  ini  use  "Os  *«iiM 
slsrout  \ rcoisr 

I AFLl C 1 7 Jw»  <*  t A-M,C-i  I 

REAL“d  XL  :o.  1C1) 

I A TeG£R  KPLS) 


JCNNUN  /?Ne/ ARAC.SAAO.SPA.SRa.  *01! 


EP$«ii.c*n«*<Mjui» 

ARAC'C £*C 


JNE-1  .0«0 

mo-o.c«o 

8i:\:8t}« 

3 ) • LT  *0 ) 
ifikpl  i)  .£  j.  :i 

E^A-dRAO-ARAC 

HftiPi*  I .lS.l  I “*11*16.2)  KPtl) 
CRMAttn  , 10X,*  *3015*  F*<5*LS.M  CP 
CO  5 1 • 1 * N 

IlM- I- 1 
MEJM/EN'L 

X 14,1)  «ARADMiftAO-AKAC  MV 
X IS, ! 1 -V 

x (o.n  «> 

11ER«  ‘ 


tVPe  Ml.'.'l 


c :hs 


VrER«  l 
M *P<  1 ) .6  J. 

na&LE*  is  . 

C 1*  L C .0  «0*PO  I i 
CALL  1N7VH.7. 


_)  CO  TO  ICC 
r TVP  e U JNMJLA* 
iMLiAAC-AAAOJ 
tfPA  I 


0 ISx 


__  l nTv l 1 1 

i"il  y.  ll-XU  ,N  Ml  CSE-TOt  SMlSPA-SRE) 
MLLVli.C.i) 

MULVM.tt.S  J 
iNTvi«»,io,e?iAi 

NnxtI,M*SA6*XU.lM$RA 

tol 

S M x(  11, N) 


CALL  “ULVM.tt.l2) 


30 


) ITS*. A. » 


CALL  t \ Tvt  li.  U.  3MAI 
CA»X(13,N| 
o*ci“C“-o;*c3 
A.<CMC*-Ctt»C3)/0 
6»LCi«Ctt-c:*C5)/o 
CONTINUE  . 

1MKPM)*C0.1)  1* 

PCRMAT15*.  IlO.IPCeCO.J) 

CALL  *Ut$l7,M,A) 

CALL  “ULSL13.  L5 . 5 I 
CALL  AOCVl  1 A,  15.  IS 
Call  slsvi  15.  u.  la) 

S *SK8*  X L l.*«»-Xtl6.NI 
CALL  AOCSl It, L4.S) 

2A-X13.1M  MAAAU*  IcNi-POIS  ) *X  l 16 , 1 » 
CALL  *ULi< 1 1. i 7. d I 

call  suea“,i3.ARACi 
call  PUL  5Llo.i3.AI 
Call  aoovl  l7.ii.iii 
]«-! CNE • c C f S A 

n, 


!all  Su6vl  i j.  j.  in 


PCCCCMO 

SSSSS<"S 

Looccs  ;o 

RCCCCS  K 

rc  " 


SCOCCtCC 
ROCCOft  10 
PC?  06  .3 

mim 

PCCCCcSC 
R0  CCOttO 

imtu 

RCCCOe^O 
RCCCC7CC 
RCOC  C 7 1 C 
RCCC07. j 
PtfoCOT  -J 


30783 

SjjlcS 

!Ce  1C 
:cs:c 
.a*5oe  :j 

SC0C0d-0 

RCCCCMC 

rocCcea: 
RCCCCt  7$ 

roccces' 

sccccesc 

RGCCCSCC 

rccccsi: 

RCC?CSiO 
RCOCCS -C 
ROD  COS  *0 
RCOCOSSO 

RO CcCS'C 
RCD  ( 1000 

SS8HSi4 

ROo£l£*c 

RCCCICfC 

PC0C104J 

m\u$ 

RCCC1C‘C 

ROCC11CO 

ip? 

R0CC1MC 
RO  CC  1 1 SC 
RCCCllaU 
RODCll/0 

Rcccueo 

mm: 

RCCC  121v 
pcdck; 
rcoci;: 
Retell  “ 
rccci:  : 
ro  cc  i;« 

RCCC1C7 
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Jftk.iVKil7-1#’ 

ifUPi  3 ).lt.g»  a 


I f<KP( 3 l.LT.G  ) GO  TC  215 
IFIITSR.CT.KP  li)  ) GJ  TO  200 
*0  CONTINUE 

CALL  0UFV(4,  14) 

If  (KP(  1)  .£0.2)  X(1<,,1)*CNE 
CALL  DI VV  ( 17,14,1  «) 

CALL  suesue,  16.  A) 

Call  ons na.c.a) 

If (KP( 1 ) .E0.2)  GO  TO  15C 
GC  TO  20 
100  A-ZSRC 

C»(0N6-P0!S  J/X( 1,1) 

CALL  MULV(1,2,7) 

CALL  NULV(4,7,a> 

CALL  I NT  V ( 3 , 10 , 2MA ) 

|lM*X(  10  ,N  )*SRo«( *11 .N >-X(l, 1) ) 

5 0**0*  $UM*x(2,i)-X(3»N) 

150  CALL  1 NTV(fc,ll,8MA) 

0EN-8RJC* (OAE*POiS)*Xlll,N) 

CALL  MUL  V (1,6,11) 

CALL  INTVdl,  13,  3HA  I 
0EN-0EN*C«X (13, N) 

8-SUM/DEN 

CALL  1 N T V ( 1 ,7 > 8MA  ) 

CALL  INTV(6 ,11, 8MA) 

GC  TO  30 
215  93*62 
3 2*6  1 
31*6 

2UH*0AeS( e 1-92 1 *0  4 8$ ( 62- 63 ) ♦ CA8S< 8 3-9  1) 
01V-OA3S(ai)»0AeS(32)*CA3S(33 ) 

CRIT  *ZUM/C IV 

1 P(CRI T.LT.SPS)  GC  TO  200 
GO  TO  40 

2 CO  CALL  ACC  V( 17 ,16,194 

If (wP( 3)  .LT.O.AND ,KP(4 ).EQ.O)  k.RITE16,201)  ITER.EFS 
201  FORMAT  <//,20x,23.  ' ITERATIONS  R6CUIRE0  WITH  EPSJLCN  » *,1PE8,1) 


If  ( KP(4  ) .i£  .1 ) »R IT S (o  , 2u4 ) 

204  FORMAT)//) 

IF (KP(4).LE.1>  HR  IT  £( 6 , 2C  5 ) 

205  FORMAT t 23X. 'RADIUS' , 10>, 'THICKNESS* . 5X. • GAMNA  OFEGA  SQ*,7X, 
1 • EE  ALPHA  TEE*  ,7X  ,*  SIGMA  RAO IAL ' , 6X , • SI GM A CIRCUMF'i 

N$KIP«KF(2) 

00  210  I »I,N,NSXIP 
J * 17  NS  K I F 

IF<KP<4)  ,l£.  1)  HR  IT  E (6,  211  )J,  X(  4,  I)  ,X(  1 ,1)  ,X(  2,1  ) , 

1X( 3,1 > . X(16  1 J , X(  19, I ) 

210  CONTINUE 

IF(KPU).GT.l)  KP  ( 5 ) * IT  E R 

211  FORMAT  ( 110, 1Pa£  19.5) 

RETURN 

ENO 

THIS  IS  THE  START  OF  THE  ANCILL  ARIES 
SUBROUTINE  40CV ( N 1 , N2 • N3 I 
REAL*  8 X(2C. 101). S 
INTEGER  KP(5) 

COMMON  X ,N, KP 
00  1 1*1, N 

1 X (N3, I J-X(N1» I) *X(N2,  I I 
RETURN 
ENC 

SL6R0UTINE  SU3VIN1  ,N2  ,N * ) 

REAL* 8 >(20, 101), S 
INTEGER  KP  < 5 > 

COMMON  X.N.KP 
00  I I * 1 <N 

1 X (N3,  I I-X(N1, I)  -X  In2»  I ) 

RETURN 

END 

SueROUTINE  MULV ( N 1 ,N2 ,N  3 ) 

REAL*  3 ><20, 101). i 
INTEGER  KP ( 5 I 


ROCC 1260 
RCQ  C 1 2 * C 
ROCC  12C0 
RCDC13  10 
R00C132C 
ROC  C 1 3 30 
ROOC124C 
ROCC1360 
RCCC13eJ 
ROC  C 1 3 70 
RC001360 
ROCC 1350 
RCDC14CC 
RODO 1 4 10 
P COO  14  20 
RC0C1F2C 
PCCC 1440 
RC0C145C 
R0CC146C 
RCCC147C 
R00014  6 C 
R0001450 
RCCC15CC 
R0DC1S  1C 
PCC01520 
ROD  C 1 £ 20 
RC0C154C 
ROCC  1 5 SO 
RCOCIStC 
ROCO  1£70 
FC0C1530 
ROC  Cl  SIC 
R00016C0 
RCCC161C 
ROD  C 1 6 2 C 
PC CC163C 
RCO  C 16  40 
PCCClcSC 
ROC  C It  6 0 
RCCClt  7C 
pccciteo 

RCDClfcSO 
ROC  C 1 7CC 
FOOO  17  10 
RCCC1 72C 
ROC  C 17  2C 
RCC  C174C 
ROC  Cl  7 f fl 
ROCC176C 
ROC  C177C 
RCDC176C 
R00C1750 

pcocieco 
RCDO  61C 
ROC  C 16  20 
RCDC162C 
ROC01&*0 
RCDCl&SO 
RCC  C 1 6 6 C 
R0DC16  70 
RCCC1  6 EC 
R0DC1C5Q 
RCDC1 SCO 
PCDC15  10 
RCDC152C 

fccc  isk 

RC0C1S4C 
ROCC  19!0 
RCC  C 15 60 
ROD  ( 1 SIC 
ROD  C 15  $0 
ROC  C 1 5 eC 
ROOC2CCO 
RC002010 
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THIS  PAOS  IS  BEST  QUALITY  FRACIlCAttU 
ERUM  CXtt'Y  iUvMStiEL  10  UUC  — 


CON MCN  X.N.KP 
CC  1 l *1.N 

l X<NJ  .1  I o<M,1 1«X<N2  , 1 1 
RETURN 
I NO 

SUBROUTINE  CJvv(.NltN2,M3) 
real«b  xico.ion.s 
INTEGER  « i=  * 5 I 
CCNHG.N  x,.\.*P 

co  i :-i.n 

I X(N3.  W»X!M,II/XIN2,I> 

RETURN 

!iif5ra!=.'?sfui-NiiS‘ 

INTEGER  XP  ( i) 

C ::**•  ( x.n.kP 

O'C  V * i »N 
I XI 
RfY 

1u?R0UT1N=  <U3S(N1.M2 .SI 

R?AL"«  X» 30,1011. i 
INTEGER  «Pt5l 
COMMON  x.n.kP 
00  1 !«1,N 
i x^N2^:>«xiM,n-s 


NO 


ROUTINE  ',UlSIMv.N2, 
IAL  «B  X!20.l3»J.a 
INTEGER  N?  ( il 
COMMON  »,N.is.P 

00  1 ! • 1 iM 

X!N2.l  » «MN1  ,1)  •$ 
TlTURN 


NO 


S ) 


ROUTINE  0Iy«<Nl.N2.S) 
L«F  XC0.10U.5 
INTEGER  .P<5) 

CCMMON  X,N,kP 

oo  i :»i.n  . , 

x(n2, ; » •*( m ,i i/s 

RETURN 

END 

SUBROUTINE  RRIVINI.I.JJ 
REAMB  Xt’T.lJD.i 
INTEGER  *P<  D 
CNMCN  ),,N.*P 


IP  5ECQNC  JRCUW'NT  *..jALS  c.  00  oirectly  tc  return 
IF  SECOND  A.GL-EM  sMALS  1, PRINT  Th6  VECTOR. 

IF  S ICON  C AACU*«NT  cCJALl  i.  JN  « N T ThE  I JEM  ITY  A AO  THE  VECTOR. 

IF  SECJNO  ARGLMcNT  t»jALS  J.  PRINT  tHs  vjCTCR  NUMBER  ANC  THE  VECTOR. 

IF  $ ECONO  AROLVEM  ^vLALS  v.  PRlvT  NUMBER.  IDENTITY,  ‘ 

IF  SECOND  ARGUMENT  cGOALS  f.  PRINT  IDENTITY  CM.Y . 

IFI  I .SO. 01  00  TO  . J 

in  i.eo  .i  i oo  to  i 

$ W i 
Y>\YM'A\  8?  \°c  $ 

1 OC  *1  K-l.N 

WR  ITE I6i?!  K.XINl.K) 

FORMAT (20X,  Il.lPfcCC.S  i 
RETURN 

R?tTE(6.;iI  J 

FORMAT! //,JCX,'  VECTOR  a I Th  IOENTITY 
GO  TO  \ 

NMTe!6.31)  N 1 

FCRMAT! //,  30X,  • VECTOR  NUMBER  • • 1 3 . ' 

GC  TO  I 

A NR  IT  E < a « 4 1 ) M.  j 
41  FORMAT!  //  ,20X  .*  VECTOR  X',I2, 

GC  TC  1 


iNC  VECTOR. 


9 

10 


31 


.15.  * FOLLOWS  * • I 


FOLLOWS**  ) 


RQCCCCCO 

lo^c'mc 


SL. 

RO  C C * v 
RCCCCMO 
RCOCCOeO 

poc£lcsS 

RCOCCC  ‘0 

Rcrcnco 

rcdccmo 

RCDC2KC 

vmm 

r?c5Ii!c 

ROOOC'.TJ 
R00C21E J 
RCCCCIM 

roccccc; 

cc2!:o 


tuw  v . . % -j 

mm 


WITH  IOENTITY  '.IS,*  Ft  LL  C WS  J • I 


FCDCC; 

$r 

•Cwi... 

RCCOCCtD 

RCOCCITO 

i&nui 

RCCC23CC 

m&i 

rcco::  ? » 

ROOCC j»C 
RCCCC3  ‘.C 

roc  c c r c c 

RCCC23  ;c 
PQCC:?E3 

8S«i^S 

RCOCC-iC 

RCCCICO 

mm 

RCOCC-cC 
ROCCCMC 
RCCOCML 
RCCOCVSO 
ROCC25CC 
RCCCC3 IC 
ROOCCfk 
RCC  C * ! 20 
PCCOCE^O 
RCCCCE iC 
RCCCC3 cO 
R CD 005  T C 
RCCO  0 f 60 

wmvc 

FCDOCt  10 

RoE^e  iC 

ROCCC  t S3 
ROOCCHeO 
RCCCCt  TO 
ROCCCctC 
RCCC24S0 
P00027CO 

m\\ 

RCCO 1 7 tO 
FCOC075O 


- S - 


— 


',15,'  HAS  EEEN  GENERATEC.'J 


5 WRITE (6.51)  J 

51  FCRMAT(//,30X,*  VECTCR  WITH  IOENT  ITY 
GO  TO  10 
ENO 

SUBROUTINE  CUPVIN1.N2) 

REAL*8  X 120 , 10H  ,3 
INTEGER  HP ( 5 > 

CCMMON  X.N.KP 
CO  1 I*1,N 
1 X (N2.I ) =XIM1, I) 

RETURN 

ENO 

SUBROUTINE  INTV<N1,N2  ,$) 

C KP15)»C  CAUSES  MILNE  INTEGRATION  TO  SE  USED. 

C OTHERWISE  TRAPEZQ  1CAL  INTEGRATION  IS  LSEC. 

REAL*  8 XI  20.1  Cl  ),  S.  AOC  ,N  I NO  ,.NTNO  , f I \ Q , THT  C . EN  , R , P 
INTEGER  KP(5) 

COMMON  X.N.KP 
I tlKPl  5 iJ'ie  .0)  GO  TO  10 
EN»N-1 
EN*1 .0*0/EN 
MNO»EN*S. 0*0/2.  aO*1 
NTNO  *EN* 1 .90  + 1/ 2. AC*1 
PI  V0*EN*5. 0*0/2. <*C*l 
TWT0*EN*1. 30*1/2. MC*1 
R*EN/  2.  AO*  1 
X ( N2.  1 J»0.0*0 

X <N2, 2)*N2-NC*X INI .1 1 *NTNC*X(N1.2 )-F I V0*X<N1 ,3)*X(N1 ,A )*P 

NM3*N-3 

CO  l K*1,NM2 

KP1*K*1 

KP2*K*2 

K P3*K*3 

acc*thto*(x  ini.kph *x (ni ,kp2>  j-r*<x<ni,khx<ni,kp3)  j 

1 X(N2«KP2)*X(N2«KP 1 ) * A 0 0 

XIN2iN)»X(N2.  N— 1 )*>IN0’*X(N1.N)*NTN0*X(N1»N-1I-FIV0*><N1  ,N-  2 ) 
1*)(N1 ,N-21*R 
CALL  ML  L S I N 2 , N2 • S J 
RETURN 

10  CONTINUE 

X (N2» 1 )*0.C*0 
P » 2* } N- 1 ) 

OC  11  1*2, N 
J-I-l 

11  X(N2,n*X(N2,JI*X(M,I)/F*X(Nl,J)/P 
MULSIN2.N2.SJ 


RETURN 
ENO 


RCCC27EC 
PCCC277C 
PCC027E0 
RCDC279C 
RC0C2ECC 
FCCC2S  10 
RCOC2E2C 
R0CC2820 
RC002SA0 
RCCC2E  SO 
R00C2etC 
RCCC2E7C 
R0DC2EEC 
RC  DO  2 6 90 
ROOC29C3 
RCDC29IC 
ROC  C 29  20 
RCDC292C 
RCCC2SO 
RCDC29  50 
RCC C296C 
ROC  C 29  70 
RCC (2  9 EC 
RO  C C 29  90 
RCC  C3CC0 
ROC C 20 10 
RCCC2C20 
ROCC2C2C 
RCDC2C90 
RQCC2CS0 
PCDC3CE0 
RCDC2C  TC 
ROCC2CEO 
ROC  C 20  90 
RCDC2  ICO 
RCC  C 2 1 1 C 
RCC C2  120 
RCCC212C 
RCCC2HC 
RCC  C 2 1 SC 
RCCC31EC 
RCCC217C 
PCCC2  ICC 
PC0C319C 
RCCC22CC 
PCCC2210 
RCCC222C 
RCOC222C 
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Appendix  B 
Sample  Problem 

A disk  rotating  ad  7200  rpm  and  composed  of  a medal  having  a 
a specific  weight  of  0.283  pounds  per  cubic  inch,  is  0.85  inches 
inside  diameter  and  5.15  inches  outside  diameter.  The  radial  stress 
at  the  inside  radius  is  zero  and  that  at  the  outside  radius  is  22,000 
psi.  The  thickness  varies  with  radius  according  to  the  law 

t 3 0.1493  r-0,4^  (all  dimensions  in  inches) 
and  the  temperature  change  (from  the  zero  stress  condition)  is  given  by 
T = o0  - 1.6  r^. 

Take  E 3 29 ,000,000  psi  and  a 3 5.7  io-c  /°P  and  determine  radial 

stress  (a  ) and  circumferential  stress  (a.)  as  functions  of  r. 
r 9 

This  problem  illustrates  most  of  the  capabilities  of  RODISK. 
3ecause  of  the  special  nature  of  the  thickness  variation,  i.e,  a 
power  relation,  an  analytic  solution  may  be  established  so  that 
the  accuracy  of  the  RODISK  solution  may  be  evaluated.  Results  of 
such  evaluations  are  given  in  Table  1 of  the  body  of  this  report. 

There  it  may  be  seen  that  accuracy  far  better  than  engineering 
considerations  require  or  Justify  may  be  obtained  by  taking,  say, 

N 3 101  and  KP(3)  * 25.  so  that  in  3.1  seconds  RODISK  returns  to 
the  calling  (i.e.,  input)  program  results  with  a maximum  error  of 
0.01  % or  less.  The  tabulation  'which  follows  shows  output  with  N 
3 101  and  KP(3)  3 -6,  resulting  in  27  iterations  and  taking  3.3 
seconds.  Accuracy  is  better  than  .006jS. 
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Alexandria,  VA  22314 

2.  Library  2 

Naval  Postgraduate  School 

Monterey,  CA  93940 

3.  Dean  of  Research,  Code  012  1 

Naval  Postgraduate  School 

Monterey,  CA  93940 

4.  Chairman  1 

Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  CA  93940 

5.  Professor  J.  E.  Brock,  Code  69Bc  15 

Naval  Postgraduate  School 

Monterey,  CA  93940 

6.  LCDR  Robert  E.  Brown,  USN  2 

Naval  Amphibious  School 

Coronado 

San  Diego,  CA  92155 

7.  Mr.  Charles  Miller  1 

NAVSEA  Code  0331 

Naval  Ship  Systems  Command 
Washington,  DC  20362 

8.  Mr.  Richard  Carleton  1 

NAVSEC  Code  6146 

Naval  Ship  Engineering  Center 
Washington,  DC  20362 


